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We study quasi-static perturbations in a cosmological background in the Dvali-Gabadadze-Porrati 
(DGP) braneworld model. We identify the Vainshtein radius at which the non-linear interactions 
of the brane bending mode become important in a cosmological background. The Vainshtein ra- 
dius in the early universe is much smaller than the one in the Minkowski background, but in a 
self-accelerating universe it is the same as the Minkowski background. Our result shows that the 
perturbative approach is applicable beyond the Vainshtein radius for weak gravity by taking into ac- 
count the second order effects of the brane bending mode. The linearised cosmological perturbations 
are shown to be smoothly matched to the solutions inside the Vainshtein radius. We emphasize the 
importance of imposing a regularity condition in the bulk by solving the 5D perturbations and we 
highlight the problem of ad hoc assumptions on the bulk gravity that lead to different conclusions. 
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INTRODUCTION 



The acceleration of the late-time universe is one of the most important problems in cosmology. Within the framework 
of general relativity, the acceleration is supposed to be caused by unknown dark energy. The simplest option for 
dark energy is vacuum energy, but it is hard to explain why the vacuum energy is so small compared with the 
J> . prediction of particle physics. An alternative to dark energy is provided by models where large-distance modifications 
0^ ' of gravity explain the acceleration. Probably the most widely studied example of a modified gravity model is the 
, Dvali-Gabadadzc-Porrati (DGP) brane- world model in which gravity leaks off the 4D brane into the 5D bulk spacetime 
' The 5D action describing the DGP model is given by 
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!>■ ; ^^TT^ d^x^'-^>R+ / d^x^R+ / d^xV^-Cm + — / d^x^K, (1.1) 

, J ^^^4 J J 1^ J 

where Cm is the Lagrangian for matter on the brane, K^^i, is the extrinsic curvature and K — K^^. The transition 
from 4D gravity to 5D gravity is governed by a crossover scale r^. 
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, . - which is the only parameter in this model. A striking feature of this model is the existence of a solution where 
rS [ the acceleration of universe is caused entirely by gravity without introducing the cosmological constant 0. In this 
• solution the Hubble parameter approaches a constant, H '^/I'c, at late times, mimicking the cosmological constant. 
■ " " ' This self-accelerating solution has attracted significant interest recently 0]. 

Unfortunately, it has been shown that the self-accelerating universe contains a ghost 0, [^, 0, 0, H, Q • The existence 
of the ghost was shown rigorously on a de Sitter spacetime by studying linearised gravity. Recently, however, there 
arc some claims that the non-linear interactions obscure the conclusion on the existence of the ghost p^ . [ll| . It 
has been recognized that the non-liner interactions of gravity in this model arc much more subtle than 4D general 
relativity [1, [5[[l3,[H, [13, E^, [lBi[l3, SiIUjH^I- The reason is that the graviton contains a scalar degree of freedom 
and the non-linear interaction of this mode becomes important on much larger scales than the usual graviton. This 
is analogous to the massive gravity model, where a helicity-0 mode becomes strongly coupled on very large scales 
for small graviton mass [2l[ . In the DGP model, the scalar mode is a mix of the helicity-0 mode of the spin-2 5D 
graviton and the spin-0 mode called the radion p, 01 • Physically, the scalar mode describes the bending of the brane 
in the bulk [J, [^, Il8| . It was shown that the non-linear interaction of the brane bending mode becomes important at 
the so-called Vainshtein radius r^, = ("Tgr^)^/^ where fg is the Schwarzschild radius of the source [l3|- If we want to 
explain the late-time acceleration, we should require Tc ~ . 

One argument against the validity of the linearised analysis is that, for cosmology, is roughly the Hubble scale 
today Tg ~ , then the Vainshtein radius is also the horizon scale , which may indicate that the linearised 

cosmological perturbations are not valid p^ . However, most of the literature so far studied perturbations around 
Minkowski spacetime. It is still unclear what is the Vainshtein radius in a cosmological background. This is an 
important question to be addressed because the ghost exists in the self-accclcrating solution where the Minkowski 
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spacctime is not a solution. There is no ghost in a Minkowski brane in the DGP model. Thus it is important to study 
non-linear interactions of the brane bending mode in a Friedmann background. 

There is also a claim that the perturbative approach cannot be applied in the DGP model [IH. This argument is 
based on the Schwarzschild solution obtained in Ref . [2^ , which does not recover the linearised solution in the region 
r > r». However, this solution is obtained by closing the equations on the brane by imposing ad hoc assumptions on 
the bulk gravity. In Ref. p3l |. it was shown that it is crucial to impose a proper boundary condition in the bulk to 
determine the behaviour of gravity on the brane. 

The aim of this paper is to study the non-linear interactions of the brane bending mode on a Friedmann background. 
We build on Ref. |23| where linearised perturbations are solved properly in 5D spacetime. We extend the analysis 
of Ref. [2^ by taking into account the second order effects of the brane bending mode. Then we study whether 
the linearized cosmological perturbations can be smoothly matched to the solutions inside the Vainshtein radius. It 
should be noted that non-linear interactions on a Friedmann background were studied in Ref. [l6| assuming spherical 
symmetry and the modified Vainshtein radius was identified. We will confirm their result by properly solvin g th e 5D 
metric perturbations without closing the equations on the brane in an heuristic way in the same spirit as Ref. [23] . For 
this purpose we closely follow the approach of Ref. [l9l |. which studied weak gravity on the Minkowski background. 



II. QUASI-STATIC PERTURBATIONS 



A. Basic equations 



In this paper, we focus on weak gravity sourccd by quasi-static matter fluctuations in a cosmological background. 
This analysis can be applied to describe the metric sufficiently far from a local source located in a cosmological 
background. We can also study the cosmological perturbations on sub-horizon scales in the matter-dominated era, 
which are relevant for the structure formation problem. 

The first order metric in the bulk is given in a 5D longitudinal gauge by 

ds^ = -(1 + 2A)N{t,y)^dt^ + (1 + 2n)A{t,yf5^jdx'dx^ + (I + 2Ayy)dy^, (2.1) 

where 

A{y,t)^a{t)ilTHy), Niy,t) = 1 T H (^l + y, (2.2) 

are the solutions for the background metric 01- Note that the (y, i)-component of the metric can be neglected for a 
static source. The Hubble parameter is determined by the Friedmann equation and the continuity equation: 

^E=H^^^p, p + m{p + p)p = Q. (2.3) 

The solution with — sign is called the normal branch solution and the solution with + sign is called the self-accelerating 
solution because there is a de Sitter solution even without any kind of matter on the brane 0. The self-accelerating 
solution attracted significant interest as a model for dark energy from the large- distance modification of gravity. 

In the 5D Longitudinal gauge, the brane is not located at y = [l^]. Then it is more convenient to move to a gauge 
where the brane is located at y = 0. We perform a gauge transformation y y — r^S,, where ^ is a scalar function 
describing the perturbation of the brane location, which is often called the brane bending mode. The resultant metric 
has the form 

ds^ = -(1 + 2^)N{t, yfdt^ + (1 + 2$)A(<, yfSijdx'dx^ + 2r^ip.,dx^dy + (1 + 2T)dy'^ . (2.4) 

At first order, is identified as the brane bending mode ^. We are interested in perturbations well inside the 
horizon. Thus we will neglect all subleading terms suppressed by aH/k <C 1, where k is the 3D wavenumber of the 
perturbations. Within quasi-static approximations, time-derivative terms can be neglected. We also neglect terms 
like {A' /A)^' where prime denotes a derivative with respect to y. This is based on an assumption that $' ~ fc$. 
This assumption will be justified later. Although we are dealing with the linearised metric perturbations, it has been 
recognized that second order terms of ip can be important on larger scales compared with the other second order 
contributions [1, [^, [H, [l^. Thus we only keep the second order terms for ip. This assumption will also be verified 
later. 
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Under these assumptions, the 5D Einstein equations are given by: 

^ 1 + 1 -^(2^ + ^) W + 4 [(^V)^ - (V.V,^)^] = 0, (2.5) 



y A"^ \ A N J ^ 2A'^ 

s(^)Gy = + 2$'), - ^ [(vv)(v,v,^') - (v.^)(vV')] = 0, (2.6) 

Si^)Gl ^ 3$" + -lv2$ + ^(r - r,^') - 2^ W + 2^ [(VV)^ - (V,V,^)^] = 0, (2.7) 

1 r / A' N'\ 

<5(5)G} = -7^(V*V, - <5]V2)($ + + r - r,(^') + 5]{^" + 2$") + ^(V^V^ - 5]V^) ^ + ^ W (2.8) 



-% [(VV)(V'V,^) - (V, VV)(V'Vfe^)] + U]^^ [{V\f - [VkVivf] = 0. (2.9) 

For the spatial components S^^^^Gj, the trace of the equation gives 

9 / A' N'\ 1 

-^V2($ + VI/ + r - r,/) - ^ (^^ + ^ j VV + 3(*" + 2$") + [(VV)' - (V,V,^)2] = 0. (2.10) 
On the other hand, taking the divergence of the traceless part of J^^-'G* , we get 

5($ + VP + r - r.^') - ^ (I + f ) VV + [(VV)^ - (V,V,^)^] = 0. (2.11) 

The existence of the brane imposes the junction condition at the brane, that relates the extrinsic curvature with 
the energy-momentum tensor on the brane 

^2 

A'^.-if5M---yT^. + rcG^.- (2.12) 

We should note that due to the induced gravity term, the Einstein tensor appears in the junction condition. The {t, t) 
component of the junction condition Eq. (|2.12p gives 

Av2$ = _K^5/3+^VV--*'- (2.13) 



The spatial components give 



$ + * = ip, (2.14) 
*' + 2<I'' = 0. (2.15) 



Solutions in the bulk 



Let us first solve the perturbations in the bulk. Combining Eq. (|2.5p and Eq. (|2.1ip . <I> and F — Vcf' are written in 
terms of \l/ and ip: 



1V2^ r, f^A' N'' 

<I> = v& -I — 2 1 

^2 2 ^2 ^ 2A'^\ A N J ' AA* 

V2 , 1V2 

-^(F-r.v.) = -- — V 







2A2 








2A2 





(2.16) 
(2.17) 



Consistency between Eqs. (|2.6p and (|2.16p requires 

(y9' = 0, *' + 2$' = 0. (2.18) 



The latter is consistent with the junction condition Eq. ()2.15p . Then substituting Eqs (j2.16p and (|2.17p into Eq. (|2.7I 
and using Eq. ()2.18p . we get a wave equation for vf 

V2 ^ /TV' 



*" + ^*H^^VV = 0. (2.19) 
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By performing a Fourier transformation, the solution is given by 

N' 



^ = * - —rc^ = 



ciilTHy)^'''''" + C2{lTHy)^^/''" , (2.20) 



for a given fc, with our approximation k/aH 3> 1. We impose the regularity condition in the bulk so that the 
perturbations do not diverge at y ^ cxd in the self-accelerating branch, and y = l/H in the normal branch. This 
means that we take C2 = [23|. We should note that the regularity condition verifies our assumption that the terms 
like {A' /A)^' can be neglected compared with the terms like W'^^/A^, with our approximation k/aH ^ 1. 

C. Equations on the brane 

Now we impose the junction conditions. From Eqs. (j2.20p and (|2.16p . it is possible to show that 

£!^A<i,«:^$, (2.21) 

Tc arc a 

for perturbations whose physical wavelengths are shorter than Tc, krd a 3> 1. Thus we can neglect $' in the junction 
condition Eq. (|2.13p . Then the projection of Eq. (|2.5p on the brane and the junction conditions Eqs. (|2.13p and (|2.15p . 
provide a closed set of equations on the brane for <f>, ^' and Lp. The effective Einstein equations are written as 

4^V2$ = -K^^p+^VV, (2.22) 

+ $ = (p, (2.23) 

and the equation of motion for tp is given by 

mt)^V + J [(VV)' - (V,V,^)2] = nlSp^ (2.24) 

where 

*)-l-f (2^ + ^)^l±2H.. (l + J,). (2.25, 

Here the + sign corresponds to the normal branch and the — sign to the self-accelerating one. 

III. SOLUTIONS ON A BRANE 
A. Linearised solutions 

by negl 

perturbations are easily obtained as 

(3.1) 
(3.2) 

V o/ji / 

which agree with the solutions obtained in Ref. [l6|, [2^ . 

The linearised equations can be described by a Brans-Dicke (BD) theory. The perturbed Einstein equations in the 
BD theory are given by 

<5G^. = -(V^V,-ff^,V2)(^, (3.3) 
and the equation of motion for the BD scalar is 

= T^^P^ (3-4) 



lectin^ 


I the second 


orde 






1 


— * 










" 3/3 








— * 
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where u) is the BD parameter. Comparing Eqs. (|2.22p - (|2.24p with Eqs . (13.31) and (|3.4p . we find that the branc bending 
mode acts as the BD scalar and the BD parameter is given by [l^, [2^ l26j 

u^ = \{^-V). (3.5) 

The sign of /3 is directly related to the existence of the ghost in de Sitter spacetime. In the self-accelerating branch, 
/9 is negative for Hr^ > 1/2. In the BD theory, the BD scalar has the wrong sign for the kinetic term if w < —3/2, that 
is /3 < 0. The condition that f3 is negative is given by Hrc > 1/2, which is precisely the condition for the existence 
of the ghost in de Sitter spacetime, as was shown in Refs [1, [1, @|- On the other hand, in the normal branch, /? is 
positive and we expect no ghost in this branch of the solutions. 



B. Spherically symmetric solutions 

Next, we study the effect of the second-order contributions of ip. Unfortunately, it is not easy to solve the equations 
for ip with the non- linear interactions. Thus we assume spherical symmetry to simplify the problem. The equation 
for If (|2.24p is then given by 

where 

E^2rl[U'^]\r, (3.7) 



r 



dr 



in agreement with Ref. [19[ in a Minkowski spacetime. Let us consider a source localized in some compact region. 
Then it is possible to integrate the equation to get 

3/3(p + S + ^ = 0, (3.8) 
r 

where 

drr^Sp, (3.9) 

is the Schwarzschild radius of the source. Hereafter, we assume rg — const, for simplicity. Taking the r derivative of 
Eq. (j3.8p gives an algebraic equation for dip /dr. Then we get a solution for dp/ dr as 

where 



2^ \ 1/3 



8r^ 

9/32 



(3.11) 



which is the Vainshtein radius for a source in a cosmological background. This is in agreement with the result of 
Ref. [l6l | , but we arrive at this result by solving the 5D bulk metric and imposing the regularity condition in the bulk, 
without closing the equations on the brane in an heuristic way. The solutions for the metric perturbations can be 
obtained as 

*^| + f, (3.12) 

On scales larger than the Vainshtein radius r > r*, the solutions are given by 
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which agree with the hncariscd solutions Eqs. (|3.ip and (|3.2p . This shows that the Hnearised solutions do make sense 
as long as we are considering scales larger than the Vainshtein radius. 

On scales smaller than the Vainshtein radius, r < r*, the solutions for and $ arc obtained as 

(3.16) 
(3.17) 

In this region, the corrections to the solution in 4D general relativity arc suppressed for r < so that Einstein gravity 
is recovered. From Eq. ()2.24|) . we can see that S dominates over the linear term in this region. This indicates that 
once (p becomes non-linear, the solutions for the metric approach those in 4D general relativity. We should note that 
(3 is negative in the self-accclcrating solution while (3 is positive in the normal branch solution. Then the corrections 
to 4D general relativity solutions have opposite signs in these solutions, as was first pointed out in Ref [lH|. By a 
simple coordinate transformation, we can check that our solutions agree with the results of Ref. [lB| . 

Even on scales smaller than the Vainshtein radius r < r^, the induced metric perturbations are small as long as we 
consider scales larger than the Schwartzschild radius r > Tg. This justifies our assumption of neglecting all second 
order contributions other than the second order terms of ip. It should be also emphasized that the (?/, r) component 
of the metric, rc^p^r, is evaluated as 

{r,^,rf ^ (7)' (7)' ^>^*' (3-18) 

^2 C'g 



{rc'P,rY - [f), for r<n. (3.19) 

The higher order terms of ip in the Einstein equations have higher order powers of rc^p^r- Thus they arc suppressed 
for r > Tg. Then we only need to keep the second order terms which can be comparable to the linear terms as is seen 
from Eqs. ([XTT]) and (jgl^ . 



C. Cosmological perturbations 

Finally, we consider cosmological perturbations in a matter-dominated universe. We define an over-density of dark 
matter as 

5=^. (3.20) 
P 

The continuity equation and the Euler equation are the same as in 4D general relativity; 

^ + iv(l + <5)v = 0, (3.21) 
at a 

^ + i(v.V)-v + iJv = -iv*, (3.22) 
at a a 

where v is the velocity perturbation of dark matter. Here we introduce time-derivative terms. In order to ensure our 
quasi-static approximation, the time-dependence of the over-density 5 should be weak, dtS <^ k6, which is indeed valid 
for dust matter. Combining Eqs. ([H?^ and ([X^ with Eqs. (P?^ and we can describe the evolution 

of the dark matter over-density. 

The non-linear terms in the equation for ip, Eq. (j2.24p . become dominant when 



Using the linear term in Eq. p.24p . (p is estimated as 



rr2 2 

.-^4 (3.24) 
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where we used k|p ^ H^. Then the condition that the non-Hnear terms become important is given in terms of S by 

P^Hr,)-'' ~ 0(1) < 6. (3.25) 
If non-Unear terms become dominant, tp is estimated as 

Then in the Poisson equation Eq. (|2.22p . the contribution of ip can be neglected and 4D general relativity is recovered. 
Thus, from these rough estimations, we expect to recover 4D general relativity for non- linear over-density (5 3> 1. 
This also means that for linear over-density S I, the second order terms of ip can be neglected and the linearised 
cosmological perturbations do perfectly make sense as opposed to the claim made in Ref. (lo| . In order to verify 
these estimations, one should solve the non-linear equations for S and ip, which is difficult even in conventional 4D 
general relativity. One approach is to consider the spherically symmetric collapse of the over-density. This was done 
in Ref. [lB|, and it was demonstrated that once the over-density exceeds 0(1), 4D general relativity is recovered. This 
confirms our estimations. 

IV. EFFECTIVE THEORY ON THE BRANE 
A. Effective theory for 

In the previous section, we find that the brane bending mode tp plays a crucial role in the DGP model. It is possible 
to understand the role of the brane bending mode in a covariant way as was sliown in Refs. [1, Q. We begin with the 
definition of the extrinsic curvature: 

i^,..= 2^(5,5^. -V^iV.-V.iV,,), (4.1) 

where ^^j^ is the induced metric, Nfj_ is a shift function and Af — \J Qyy — N ^N^^ is a lapse function. Let us first 
consider perturbations around Minkowski spacetime 

Otj^iy ^ Vt^i' + Sg^iu, (4.2) 

The lapse function is given in terms of the brane bending mode tp by N^^ ~ re's/ ^(f. Tlien the extrinsic curvature is 
given in term of by 

5K^^ ^ -r,V^V^(/j. (4.3) 

An important result obtained by solving the 5D perturbations is that we can neglect the y derivative of the induced 
metric in the junction condition because 

^yg^,u < rcV'^Sg^^. (4.4) 

Then the junction condition becomes 

(5G^^ = nl5T^^ - (V^V^ - 5m^V2)(^. (4.5) 

On the other hand, the Gauss equation in the bulk, that is the (y, y) component of the 5D Einstein equations, gives 

R- + K^^K'''' ^0. (4.6) 

Then combining Eqs. ()4.5|) and (j4.6|) . we get the equation for ip as 

3VV + rl [(VV)' - (VpV,¥>)2] = -4t, (4.7) 

which reproduces Eq. (|2.24p for static perturbations. We should emphasize that the non-linear terms for (p come from 
the non-linear terms of X^i/. Even if we are dealing with weak gravity where the induced curvature is small, this does 
not necessarily mean that the non-linearity of iC^j^ can be neglected. We should also note that the higher order terms 
in (p comes from N^N^^ which is given by {vcip^rY ^ spherically symmetric spacetime. We have shown that these 
contributions are suppressed as long as r > r^. 
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In a cosmological background, the extrinsic curvature has contributions from the background 

Ku 



N' A' 



A 



This gives an additional first order contribution in the Gauss equation [2^ 



'AT' A'^ 



(4.8) 



(4.9) 



This modifies the coefficient of the Hnear kinetic term for to 3/?. In de Sitter spacetime, this is exactly the origin 
of the ghost in the self-accelerating solution. 

The equation of motion for (p can be derived from the action 



/ 



3/?(Vv3)2 +r2(V¥>)2vV 



(4.10) 



assuming static perturbations. Defining a new field tt as tt M/^(p^ where K4 = 1/Af|, the action can be rewritten as 



S cx 



3/3(V7r)2 + _(V7r)2v27r 



(4.11) 



where A = (Af4/r^)-'^/^. In de Sitter spacetime, this agrees with the boundary effective action for the brane bending 
mode derived in Ref. [1, 0| . Thus our solution is consistent with the effective theory for the brane bending mode of 
Refs. 0,11. 



B. Effective equation on the brane 



It is also possible to construct an effective theory for ip using the effective equations on the brane. Projecting the 
5D Einstein equations on the brane, the effective equations are given by p7j 



where 



n 



24' 



(4.12) 



(4.13) 



(4.14) 



and is the projection of electric part of the 5D Weyl tensor. For fluctuations around the vacuum Minkowski 
spacetime, Gf^, is written solely in terms of ip from E g. (14. 5[) . Thus the effective equations are written in terms of ip 
except for E^^y. The resultant effective equations are (28| 



-(V^V,-g^,V2)^ = -^[5^,{(VV)'-(VaV;3<^)'}-2{(VV)(V^V,^)- (V^V„(^)(V,V"(/p)}J-£;^,. (4.15) 



Taking the trace of this equation gives 



3VV + rl [(VV)' - (VpV,¥>)2] = 0, 



(4.16) 



because E^^ is traceless. This reproduces the equation of motion for ip, Eq. (|2.24p . On the other hand, the (i,t) 
component gives 



(4.17) 



If we neglected Eu, this equation would contradict Eq. (|4.16p as is pointed out by Ref. [2S 
However, we should not neglect E^^. From the 5D metric. Eft is calculated as 



Et, 



(4.18) 
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It should be emphasized that Ett contains the second derivative of the metric with respect to y. Therefore, unless 
we solve the 5D perturbations, it is impossible to evaluate this term on the brane as the junction condition on the 
brane only determines the first derivatives. Using the solutions Eqs. (|2.16p and (|2.17p and the equations on the brane 
Eqs. (12221), (p:^ and ([221, Eu is evaluated as 

Ett = I [(VV)' - (V^V.^)2] . (4.19) 

Then it turns out that the {t,t) component is fully consistent with Eq. (|4.16p . Thus the effective equations (|4.12p 
are consistent with our solutions. This is in fact trivial as the effective equations are nothing but the projection of 
5D Einstein equations. Thus as long as we solve the 5D equations, the solutions should trivially satisfy the effective 
equations. 



C. Condition on iS^^ 

At linearised level, it was shown that the regularity condition for bulk perturbations gives a condition on E^j^^ 
which cannot be determined by equations on the brane [2^ . Here we check that this condition is not modified by the 
inclusion of non-linear interactions of (p. First let us parameterize Efj_i, as 



SE* 



1 



(4.20) 



where Stt'^e i = V'VjcJTr^; — (l/3)(5j-V^7r£;. In a cosmological background, these quantities are given by 



4SpE = 3$" + — (r-r,(p'), 
kIStte = r - Vcip'. 

Then using Eqs. (|2.16p . (|2.17p and (|2.19p . it is possible to show that these satisfy the condition 

5pE + 2^5nE = 0. 



(4.21) 
(4.22) 

(4.23) 



Note that we have already used the regularity condition to assume $' ~ and neglect terms suppressed by aH/k <C 1. 
This is exactly the condition obtained in Ref. [2^. As in the Minkowski case, we can evaluate SpE as 



i^lSpE = 



1 ± IHrr 1 



1 ± 2Hr, 



[(vV)' - (v.v,^)2 



1 ± 2Hrr 1 



H 
2IP 



1 ± 2Hr, 



K^Sp. 



(4.24) 



This agrees with the result obtained in Ref. [2^ if we neglect the second order terms of (p. The (t, t) component of 
the effective Einstein equation gives 



SG 



A'r, 



* A a2 ■ " 2a 



(4.25) 



where we used the expressions for 11^^ in terms of the extrinsic curvature 



(4.26) 



and Eqs. (j4.3[) and (|4.8p . Using the solution for 6pE, it is possible to check that this equation reduces to Eq. (|2.22p 
using Eq. (j^TM)) . 



V. CONCLUSION AND DISCUSSIONS 



In this paper, we studied quasi-static perturbations in a cosmological background in the DGP brane world. Using 
Gaussian coordinates, we derive the solutions for weak gravity by taking into account the non-linear interactions of 
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the brane bending mode. Solving the bulk metric perturbations and imposing a regularity condition, we got a closed 
set of equations, Eqs . (|2.22p . (|2.23p and (|2.24p on the brane. At linearised level the theory is described by a BD 
theory with the BD parameter given by uj ~ 3(/?— l)/2, where /3 is given by Eq. (|2.25p . We studied the effects of 
non-liner interactions of the brane bending mode assuming spherical symmetry. We find that the Vainshtein radius 
at which non-linear interactions of the bending mode become important is given by = ry/P"^ where ry is the 
Vainshtein radius in the Minkowski background. In the early universe, 3> 1, so the Vainshtein radius is very small. 
Note that in this limit, we recover 4D general relativity even at linearised level, as the BD parameter becomes large 
[25I. [26t. On the other hand, in the self-accelerating universe, Htc = 1, = 1, so the Vainshtein radius is the same 
as in the Minkowski background. On scales smaller than the Vainshtein radius, r < r, , the solution approaches 4D 
general relativity. Our solutions agree with the results of Ref. [H, [lB| in the Friedmann background, and the results 
of Refs. [13, m, 111 in the Minkowski background. 

Our equations can be applied to cosmological perturbations on subhorizon scales in the matter-dominated era. 
Although the non-linear equations are difficult to solve in this case, we can estimate the scale at which the non- 
linear interactions of the brane bending mode become important. We found that once the dark matter over-density 
becomes non-linear, the non-linear terms of the bending mode also become important and the behaviour of metric 
perturbations approach to 4D general relativity. This result is in accord with the finding in Ref. where a spherical 
symmetric collapse is studied in the self-accelerating background. Our result indicates that the linearised cosmological 
perturbations analysis does make sense in the same way as in the conventional 4D cosmology as opposed to a claim 
made in Ref. p^ . 

We checked the consistency of our solutions with the effective equations on the brane. First, we checked that our 
solutions can be derived from the boundary effective theory for the bending mode derived in Refs. [3, Following 
Ref. [2^, we also checked the consistency of our solutions with the effective equations on the brane derived by a 
projection of 5D Einstein equations. A key quantity is the electric part of the bulk Weyl tensor projected onto the brane 
E^^. If we neglected this Weyl contribution, the effective equations were inconsistent. Using the solutions in the bulk, 
we can evaluate E^j_i, on the brane. We have shown that once the contribution from Ef^i, is properly taken into account, 
the effective equations are fully consistent. Our analysis is consistent with the boundary effective action Eq. (|4.1ip at 
least for static perturbations. It was pointed out that this effective action manifests superluminal propagation if we 
consider time-dependent fluctuations around a spherically symmetric solution (29j . It would be important to extend 
our 5D analysis to include time-dependent perturbations to check the validity of the boundary effective action with 
time-dependent perturbations, and understand the causality of the propagation in the 5D spacetime. 



Our conclusion is different from that of Ref. [11|. In Ref. [11[, it is argued that the linearised perturbations, which 
by themselves are valid at r > ?% are not guaranteed to match to the solution inside r < Th,. This argument is 
based on the Schwarzschild solution obtained in Ref. 1221 . which does not recover the linearised solution in the region 
r > . However, the Schwarzschild solution in Ref. |22| is not derived by solving the bulk metric and imposing a 
proper boundary condition in the bulk. Instead they imposed a specific form of the metric and closed the equations 
on the brane. This is in fact the same as imposing an ad hoc condition on E^,j. As we have shown in this paper, the 
condition on E^i, has to be determined by solving the bulk metric and imposing an appropriate boundary condition 
in the bulk. For weak gravity that is valid for r > rg, the regularity condition in the bulk uniquely determines a 
condition for -E^j^, Eq. (|4.23p . The Schwarzschild solution found in Ref. (2^ does not satisfy this condition in the 
weak gravity region. Thus their solution is unlikely to describe weak gravity sourced by a physical local source on 
a brane. On the other hand, it is still an open question what is a proper condition on E^i, for strong gravity. An 
outstanding open question is to find a fully non-linear spherically symmetric solution that properly reproduces the 
solutions Eqs. (|3.10p and (|3.13p for weak gravity. We will come back to this issue in a separate publication. 

Finally, we comment on the ghost problem. Our analysis shows that the linearised analysis does make sense as 
long as we consider scales beyond the Vainshtein radius for a local source. Then on scales r > r» we find a ghost 
in the self-accelerating universe. Usually, we expect an instant instability of the spacetime in the presence of the 
ghost. Then the self-accelerating universe would not be a viable background for cosmology. However, in this case, 
it is not so obvious that the ghost leads to an instant instability of the spacetime classically [30| . or even quantum 
mechanically [3l| . Furthermore, non-linear interactions of the bending mode would become important if instabilities 
kick in. Further study is needed to understand the fate of this ghost. On the other hand, the normal branch solution is 
free from the ghost. Although the solution itself cannot be an alternative to dark energy, it still provides an interesting 
possibility to realize an expansion history of the universe which is equivalent to a dark energy model with an equation 
of state less than —1 [s^]. This model also provides a concrete example for the large distance modification of gravity 
Q. Our equations (|2.22p . (|2.23p and (|2.24p are the basis for the study of structure formation tests in this model. 
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